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THE EFFECTS OF CURVATURE AND VI SCOSITY ON BAROCLINIC 
INSTABILITY - A IWO-LAYER MODEL 

I. INTRODUCTION 


The near-zero gravity environment of Spacelab presents an opportunity 
to perform geophysical fluid flow model experiments in true spherical geometry. 

A dominant radial dielectric body force, which is analogous to gravity, can be 
achieved over a volume of fluid held between two concentric spheres by applying 
a large voltage across the spheres [1,2]. There is no way such experiments 
can be performed in a laboratory on the Earth^ s surface because the dielectric 
body force is weak and its effect is overwhelmed by terrestrial gravity. 

A spherical model of the general circulation of the Earth^ s atmosphere 
is being considered for Spacelab [2] . In the past, much sophisticated general 
circulation modeling has been performed in cylindrical geometries, and much 
has been learned which is relevant to the general circulation [3,4,5]. However, 
these models have serious limitations when attempts are made to extend their 
results to the spherical geometry of the atmosphere. 

Systematic scaling of the governing equations for large-scale (synoptic- 
scale) processes in the atmosphere and the laboratory models yields an identical 
set of approximate equations for both systems known as the quasi-geostrophic 
equations [6] (except that curvature terms are missing from the cylindrical 
models); therefore, the laboratory experiments are called models. These 
equations do not include external heating processes which are important on 
longer time scales, and they do not include local dissipation processes which 
are small. Using the quasi-geostrophic equations, a linear stability analysis 
of a baroclinic zonal current contained between two rigid parallel boundaries 
has been completed. Some effects of curvature and boundary dissipation were 
included. The motivation for this analysis was to guide the design of the 
Atmospheric General Circulation Experiment (AGCE) for Spacelab, This report 
describes the stability analysis and presents the results. For additional infor- 
mation on AGCE and on baroclinic instability and the two-layer model, the reader 
is referred to the references herein. 



No one has succeeded in solving analytically the equations for baroclinic 
instability on a sphere — they are not separable. To proceed, simplifications 
must be made. The analysis described in this report was performed on a 
/3-plane. The /3-plane is a tangent plane to a sphere at some latitude d on which 
rectilinear cartesian coordinates are used and on which the variation of the local 
vertical component of rotation rate with latitude is represented linearly. The 
important effect of curvature introduced in this way is known as the /3-effect. 

The effect of viscosity at the boundaries was introduced using the Ekman com- 
patibility conditions [7]. These conditions enable one to omit a detailed dis- 
cussion of the Ekman boundary layers on the horizontal surfaces while keeping 
the effect of these layers on the interior flow. No variation of the current in 
the north-south direction in the basic state was permitted. The continuous 
vertical structure was reduced to two layers within which all dependent quanti- 
ties were kept constant. This two-layer model has been used by many workers 
and has proved to be a useful theoretical model for understanding baroclinic 
instability [6,8,9,10]. 

The analysis of the effect of Ekman damping on baroclinic waves on a 
/3-plane, using a two-layer model, was first performed by Holopainen [11]. 
However, since Holopainen’ s formulation was for the atmosphere and his 
results were not presented in a form suitable for discussing the spherical model, 
it was necessary to repeat the analysis from the beginning. Holopainen dis- 
covered an unexpected destabilization of the waves by viscosity, which is dis- 
cussed in Section V. 


II. FORMULATION OF THE TOO-LAYER MODEL 


Let us consider a /3-plane on which the axes (x,y, z) correspond to 
eastward, northward and vertical directions, respectively. Let the corre- 
sponding velocity components be (u,v,w). The quasi-geostrophic set of equa- 
tions on the /3-plane is 

. ( 1 ) 


— V^p = -fV • (^ X V ) , (2) 

P - -n 


2 







where 

^ = vertical component of relative vorticity 

y = nondivergent part of the horizontal flow 

v = divergent part of the horizontal flow 
“d 

V = horizontal operator 
$ = unit vector in the z -direction 
p = pressure 
T = temperature 
p = density 

g = gravity (or some other body force) 
f = twice the vertical component of the rotation rate 
/5= df/dy. 

For the derivation of these equations and for further discussion, see References 
6 and 12. 
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since V • V =0, we can write v In terms of a stream function ib as 
— n ~n 


V = z X ( 7) 

— n — 


and 


f=z.VXv =V^jp 
“ —n 


( 8 ) 


Substituting from equations (7) and (8) Into equations (l) and (2), we obtain 


— (V^!b)+(v .V)V^i/) + j3v =-fV.v^ 

at -n ' n -d 


(9) 


-V^p=fV^i/) . (10) 

p 


Eliminating further and using equations (9) and (lO) and (3) through (6), we 
obtain 


— + (v 

at ^ ' -n 


v)v^jp + /3v 

n 


aw 

az 


( 11 ) 


at \az / '-n 




(12) 


where x = go; 7 — • 

O Z 

Equations (ll) and (12) can be solved as they are with appropriate 
boundary conditions [13] , but the solutions are In the form of confluent hyper- 
geometric functions, and computations based on such functions with complex 
arguments require much time-consuming programming. A two-layer model Is 
used. 
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In the two-layer model the fluid is divided into two layers bounded by 
surfaces, as shown in Figure 1, The surfaces at levels 0, 2, and 4 divide the 
fluid into the two layers. In each layer, the values of stream function and 
vertical velocity (and other dependent variables) are assumed constant through- 
out the layer and are set equal to their values at the midpoint of the layer. The 
midpoints are denoted by levels 1 and 3. The vorticity equation (ll) is applied 
at levels 1 and 3 with 6w/ 9z evaluated as 



(13a) 


(13b) 


where Az is the thickness of each layer. The resulting vorticity equations are 


9 f 

— + (Vj • v)v^!/)l + /3VJ =— (W 2 - Wo) , (14) 


9 f 

— + (X3 • V)V^!/)3 + /3V3 = — (W4 - W 2 ) . (15) 


Equation (12) is applied at level 2 as 


9 y /\ 2 

(i/'s - ipl) + (I 2 • V )(l/)3 - !/)i) = - ^— W2 . (16) 

^ is not a predicted quantity of the model, but is obtained by linear interpolation 
between levels 1 and 3, 


h = 2 ^ *^ 3 ) 


(17) 


5 


L 


z = 2Az = d 


4 


’*4 




3 


W2 

z-Az ‘ 


2 




1 


*0 

z - 0 0 

Figure 1, Arrangement of variables in the vertical for the 
two-layer model. 


Therefore, remembering equations (7) and (l7), equations ( 14) , (15), and (16) 
make up three equations in the three unknowns ipi, ipz, and W 2 » Wq and W 4 are 
specified by the boundary conditions. 


A. The Basic State 

We shall perturb a geostrophic baroclinic steady state specified by: 


Vi = V3 = 0 , Wq = W2 = W4 = 0 , 

(18a) 

Uj = Ul , U3 = Ug , 

(18b) 

!/>i = -Uiy , h = * 

(18c) 

U (Ug -Ul) = ^ Az , 

T 2 ^ f dz 

(I 8 d) 
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where Ui and U 3 are constant horizontal velocities in levels 1 and 3, respectively, 
is the thermal wind, and Tg(y, z) is the basic state temperature distribution. 

Equations (18) are steady solutions of equations (14), (15), and (16). 


B. The Perturbed State 

We now define perturbation quantities as follows: 


4>i= -Uiy + i/'i(x, 

y.t) , 


(19a) 

h = -Usy + i/’3(x, 

y.t) , 


(19b) 

W2 = W 2 (x,y,t) 

• 


(19c) 

The superfbc denotes a perturbation quantity. Substituting from equation (19) 
into equations (14), (15), and (16) and linearizing by neglecting products of 
perturbation quantities, we obtain 


, dx^ dyy 


( 20 ) 


a Vs , a V \ 

dx^ dy^ / 


( 21 ) 


M - 4>i) - u 


( 22 ) 


where U = 1/ 2 ( Ui + Uq) . 
m 
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Scaling shows that the direct effect of viscosity on the interior flow is 
weak, and hence viscous terms are not present in the quasi-geostrophic equa- 
tions ( l) through ( 6) . The dissipative effect of viscosity enters this problem 
through the vertical velocity forced at the boundaries by convergence and diver- 
gence of the flow. These vertical velocities are given by the Ekman compati- 
bility conditions, namely. 


at z = 0 


at z = d 



(23a) 


(23b) 


where v is the kinematic viscosity and is the rotation rate (f = 2S2). 


III. SOLUTION OF THE MODEL 


The two-layer model stability problem is now completely formulated 
by equations (20), (21), and (22) and boundary conditions of equations (23), 
We assume solutions of the form 


jj)\ = A exp i(Xx + py - \ct) , 

(24a) 

i /»3 = B exp i(?Lx + py - \ct) , 

(24b) 

W 2 = C exp i(\y + py - Xct) 

(24c) 


Substituting equations 


(24) into equations (20), (2l), and (22), we obtain 
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,,l 


ix{ (c - Ui)k2 + /3} A = ^ (C - wj) 

(25) 

ix{(c -Ug)k2 + = ^ (w; -C) 

(26) 

-i\(c - U 3 )A+ i\(c - Ui)B - C 

(27) 


where = \^ + /n^ and x = go; dTj/ dz. 

We now proceed to apply the boundary conditions of equations (23). We 

obtain 


ix{ ( c - Uj) k^ + j8 + i6) A - ^ C = 0 , 

(28) 

i\{ ( c - U 3 )k^ + /3 + i6} B + ^ C = 0 

zxz 

(29) 

XAz 

-iX(c - U 3 )A + i\(c - Ui)B - ^ C-0 

(30) 

where 6 = (fk/^ 2AAz)( t»/fi) , 



The condition for a nontrivial solution of equations (28) through (30) is 
iX{(c - Ui)k2 + j3+ i6} 0 

0 ix{(c - Ug)k2 + /3+ i6> 

-iX( c - U3) i\( c - Uj) 


J_ 

Az 

J_ 

Az 

XAz 


= 0 


(31) 
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Mialtiplying out, we have 


iX{(c - Ui)k^ + /3+ 16} |^iA,{(c - U 3 )k^ +)3+ 16} “ Ix(c-Ui) 

-(^)l-x2(c -U3){(c -U3)k2 + /3+16}] =0 



X 

, AzJ 


(32) 


Rearranging equation ( 32) to give a polynomial In c, we obtain 


k2(k2 + 2 s2)c 2 + (2(k2 + s 2)(/3 + 16) - (Ui + U3)k2(k2 + 2s^)} C + {k^iU3 


+ (/3+ 16)2 - (Ui + U3)(k2+ s2)(/3+ 16) + k2s2(Ui2 + U 32 )} = 0 , 


(33) 


where = f^/yAz^. 

Solving equation ( 33) for c gives 




(34) 


Separating equation ( 34) Into real and Imaginary parts, we have 




+ 4/326^s®] ^ cos 


V'4 i. 


(35) 
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( 36 ) 




~ ^ k2(k2 + 2 s2) "■ 


+ 4/325238] V4 gin ^ 


where 


tan I = 


2/35s^ 

(^2 _ 58) g4 + U^2 ij4(j^4 _ 4g4) 


(37) 


The condition for marginal stability is found by letting ^(c) =0. Using 
equations ( 36) and ( 37) , we obtain 


+ (k2 + s2)2s2 = -U^2k2(k2 _ 2g2) (^2 + s2)2 


(38) 


It is convenient to write equations (35), (36), and (38) in nondimensional form; 
these equations are, respectively, 


cr* 

<R(c*) 


(7* 

4(0*)=^= 


eB 


U* - 
m 




eS 




^lk*2 ^k*2 + I ^ 2Tik*2 ^k*2 + I ^ 

I ) 


M cos ^ 


(40) 


Ti ^k*2 + I j 2 Tik*2 ^k*2 + I ^ 


V4 . L 


M sin 


( 41 ) 
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II 

/ 


t ,2 2 k *4 

+ — — 
X*2 


+ ^k*2s^ = 2k*2S^l + ^k*2Sy^l - ^k*2s'' 


where a* is the angular frequency and cr* is the growth rate and where, 
K I 


(42) 


M = { B 


2 2k*4 ^-1 a\ 64 64' 

- — if + k«^k« . ^ 


+ 16 X 2048 


k*4B^ 


X*2Ta^2 g8 

m 


tan I = 


32 N/Tk*2B 



t. 2 2k*4 -1/ ' 

- — — Ta^2 


16 k*4s2 


X*2 


, *4 64 

k*4 - -yr 


m 


(43) 

(44) 


An asterisk denotes a nondimensional quantity; B = Fj tan (p/T 2 e, Ta = 

_ ^ 

eVFi'^E^, S = go AT /f^d; Fj = h/d, Fg = R/d, e = AT /AT^, E = u/fd^; and 
V n d 

h is a horizontal scale, d is the depth, R is the mean radius, and 0 is colatitude, 

AT, is the imposed horizontal temperature difference of the basic state, and 
n 

AT is the imposed vertical temperature difference of the basic state, k* and 
d 

X* are wave numbers nondimensionalized by d. U* and c* are nondimensionalized 

by fd. a* and <r*are nondimensionalized by f. 

Iv I 

IV. THE RESULTS 


Before using equations (40) , (41) , and (42) to plot the results for ranges 

of the dimensionless parameters (e, B, S, and Ta ) and the wave numbers 

m 

( X* and ju*) , we shall convert X* and p* to actual discrete zonal and meridional 
wave numbers for finite sphermal geometry. Consider our ;8-plane at latitude 
0 on a sphere of mean radius R; then. 
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9 


9 


(45a) 


2 TT 27TR cos g 
\ m 



m 

r 2 cos 6 


where m Ls the actual discrete zonal wave number. If we assume that the wave 
fills the meridional scale h, we have 



i. e. , fi* 


TT 

h 


(45b) 


We now proceed to plot the marginal stability diagrams using equation (42). 

Figures 2 through 11 show S versus Ta for several different values of m and 

m 

for fixed values of e, Fj, F 2 , and B. Each figure also shows the results for 
B = 0, The fixed parameters are varied systematically as one proceeds from 
diagram to diagram. 6 = 45® for all the figures. The Table presents a list of 
the figures and the corresponding values of the parameters. 



Figure 2. Plots of the marginal stability curves for different wave num- 
bers (m) and for e = 0. 2, F^ = 1.50, and F 2 = 2.50 [the continuous 
curves are for B = 3, 00 and the dashed curves are for zero curva- 
ture (B = 0) ] , 
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Figure 3. 'Plots of the marginal stability curves for different values of 
m and for e = 0. 3, = 1.00, and F 2 = 2. 50 (the continuous curves 

are for B = 1. 33 and the dashed curves for B = 0). 



Figure 4. Plots of the marginal stability curves for different values of m 
and for e = 0. 3, Fj = 1. 50, and F 2 = 2. 50 ( the continuous curves 
are for B =2.00 and the dashed curves for B = O). 



Figure 5. Plots of the marginal stability curves for different values of m 
and for e = 0. 3, = 2. 00, and F 2 = 2. 50 (the continuous curves 

are for B = 2. 67 and the dashed curves for B = O). 



Figure 6. Plots of the marginal stability curves for different values of m 
and for e = 0.4, Fj = 1. 50, and F 2 = 2. 50 (the continuous curves 
are for B = 1. 50 and the dashed curves for B = O) . 






Figure 9. Plots of the marginal stability curves for different values of m 
and for e = 0, 3, Fj= 5.00, and F 2 = 5.00 (the continuous curves 
are for B = 3. 33 and the dashed curves for B = O) . 



Figure 10. Plots of the marginal stability curves for different values of m 
and for e = 0. 3, Fj = 1. 50, and F 2 = 10. 0 (the continuous curves 
are for B = 0, 50 and the dashed curves for B = O) , 
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Figure 11. Plots of the marginal stability curves for different values of m 
and for e = 0. 3, Fi = 5,00, and F 2 = 10.0 (the continuous curves 
are for B = 1.67 and the dashed curves for B = O), 


TABLE. MARGINAL STABILITY CURVE PARAMETERS 


Figure No. 


0.2 


Ti 


1.50 



B 


2.50 


3.0' 


3 

4 

5 


0.3 


0.3 


0.3 


1.00 

1.50 

2.00 


2.50 

2.50 

2.50 


1.3 

2 . 0 ' 

2 . 6 ‘ 


6 


0.4 


1.50 


2.50 


1.5' 


7 

8 
9 


0.3 


0.3 

0.3 


1.50 
3. 00 
5. 00 


5. 00 
5.00 
5.00 


1 . 

2.00 

3.33 


0 

1 


0.3 

0.3 


1.50 
5. 00 


10.0 

10.0 


0. 500 
1.67 



Figures 12 and 13 show the phase speeds (c*) and growth rates (cr*) 

R I 

plotted using equations (40) and (41), respectively. Figures 12 and 13 also 

show the ordinate plotted against m for different values of S. Both plots are 

for the conditions of Figure 8 ( e = 0, 3, Fj = 3. 00, F 2 = 5, 00, B = 2. 00) and 

for Ta = 10®. 
m 


The stability diagrams can be plotted and discussed in several different 

ways. The plots chosen are those closest to the way in which the experimental 

annulus results were presented. Each figure should be considered as showing 

the results for varying AT, and S2 only, and for fixed viscosity and geometry. 

h 

Each plot assumes that e is constant. In fact, e is a dependent quantity, and its 

behavior should be determined by a complete solution of the basic state. This 

is a much more difficult problem than the one considered, but measurements 

from the annulus experiments show that e is primarily a function of the geometry 

and is only weakly dependent on AT, and . S can be considered as a measure 

h 

of the thermal driving force, Ta^ as a measure of damping, and B as a measure 

of curvature. Note that e occurs in Ta and B. 

m 

The results show that the flow is unstable within a range of S for large 

values of Ta and moderate values of B. For small values of Ta , corre- 
m m 

sponding to large dissipation, the marginal stability curves form a knee to the 

left of which the flow is stable for all S. For large values of S the flow is stable 

even for asymptotically large Ta^^. This is a short wavelength cutoff, known 

as the Eady cutoff. In this model it is due to the fact that eventually the unstable 
wavelength is too long to fit into the finite geometry. For small values of S 
below the knee there is a long wavelength cutoff. 

Note that increasing B from zero has at first almost no effect on the 
short wavelength cutoff but has a significant effect on the long wavelength cutoff. 
Even for relatively small B the effect of curvature is observed on the lower 
transition. When B > 4, the instability is completely suppressed. 


Figure 12. The phase speed versus m for different values of S 

and for the conditions of Figure 8 (e = 0, 3, Fj = 3. 00, F 2 = 5. 00, 

B = 2. 00) and for Ta = 10^. 


Figure 13. The growth rate (cr*) versus m for different values of S 

and for the conditions of Figure 8 ( e = 0. 3, F^ = 3. 00, F 2 = 5, 00, 
B = 2.00) and for Ta = 10^. 




V. SPECIAL CASES 




We have performed a stability analysis in which curvature and dissipa- 
tion are present. It is interesting to examine the following special cases: 

(l) both curvature and dissipation equal to zero (B = 0, E = O) ; (2) curvature 
equal to zero, but dissipation not equal to zero (B = 0, E # O); and (3) curvature 
not equal to zero, but dissipation equal to zero (B # 0, E = O). 

For B = 0 and E = 0 ( Ta = « ) , the marginal stability criterion equa- 
tion (42) reduces to 


S = 


8 

k*2 


(46) 


This analysis was first performed for a continuous linear vertical shear model 
by Eady [14] , who obtained the result 


k*2 


(47) 


Both results give a short wavelength cutoff only. We see that the two-layer 
model is in approximate quantitative agreement with the continuous model. 

For B = 0 and E o, the criterion equation (42) reduces to 


8 


V2 

\*2xa 

m 


8 

k*2 


(48) 


This formula was used to plot the results in Figures 2 through 11 for B = 0. 

The continuous linear shear model for the purely dissipative problem was solved 
by Barcilon [15] , and his results are in good qualitative and approximate quanti- 
tative agreement with the two-layer model results. 
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For B it 0 and E = 0, 


the criterion equation ( 42) reduces to 


= 2k*2S - ^k*2S^ 


( 49 ) 


If one performs the preceding analysis with E = 0 from the beginning, a slightly 
different result is obtained, namely. 



The stability criteria equations (49) and (50) are plotted in Figure 14 as S 
versus B for different values of m. Note the small area in which equation (49) 
envelops equation ( 50) . This means that a small amount of viscosity destabilizes 
the model in that area,. This result was first obtained and discussed by 
Holopainen [11], The finite amplitude aspects of this problem are discussed 
by Romea [16] . 

The continuous, inviscid, /3-plane model was first presented by Charney 
[13], Charney found a long wavelength cutoff, but later work [17,18] showed 
that this cutoff existed only in a limited sense. An up-to-date account of the 
continuous model results for linear and nonlinear shears is given by Geisler and 
Garcia [19] , 
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Figure 14. Plots of the inviscLd, marginal stability curves as S versus 
B for different values of m and for the conditions of Figure 8 
( Fi = 3. 00, F 2 = 5. 00, B = 2. OO) [the continuous 
curves are for vanishing viscosity, equation 
(49), and the dashed curves for zero 
viscosity, equation ( 50) ] . 
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